Majority-spin non-quasiparticle states in half-metallic ferrimagnet Mn 2 VAl 
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The density of non-quasiparticle states in the ferrimagnetic full-Heuslers Mn2VAl alloy is cal- 
culated from first principles upon appropriate inclusion of correlations. In contrast to most half- 
metallic compounds, this material displays an energy gap in the majority-spin spectrum. For this 
situation, non-quasiparticle states are located below the Fermi level, and should be detectable by 
spin-polarized photoemission. This opens a new way to study many-body effects in spintronic-related 
materials. 



Half metals display a particular type of itinerant- 
electron magnetism as well as unusual electronic prop- 
erties: they are metallic for one spin channel, and in- 
sulating or semiconducting for the opposite one P, Q. 
Electronic structure calculations based on density func- 
tional theory offer an explanation for the half-metallicity 
based on the interplay between the crystal structure, the 
valence electron count, the covalent bonding, and the 
large exchange splitting in addition to symmetry con- 
strains. The expected 100% spin polarization of half- 
metals turned out to be an excellent motivation in de- 
veloping the field of spintronics both from a theoretical 
and an experimental point of view @, S| ■ In reality many 
potential half-metallic ferromagnets exhibit a dramatic 
decrease of bulk spin polarization at temperatures well 
below their Curie temperature. In order to understand 
such a behavior from a theoretical point of view it is nec- 
essary to consider finite temperature many-body effects 
0- 

An important effect of dynamical electron correlations 
in half- metals is the existence of non-quasiparticle (NQP) 
states 0, H, 111- These states contribute significantly 
in reducing the tunneling transport in heterostructures 
containing HMF 0, H, M, E, EH , even in the presence 
of disorder. NQP states strongly influence the value 
and temperature dependence of the spin polarization in 
HMF [l@,[l3, which is of primary interest for po- 
tential applications. These states originate from spin- 
polaron processes whereby the minority spin low-energy 
electron excitations, which are forbidden for HMF in the 
single-particle picture, are possible as superpositions of 
majority-spin electron excitations and virtual magnons 
@, i, i, @. Recently we have applied the LDA+DMFT 
(local density approximation plus dynamical mean field 
theory) method (for review of this approach, see RefO) 
to describe from first principles the non-quasiparticle 
states in several half-metals 14 , 15, ljl 17, [lj|. Up to 
now, our studies were restricted to half metals with a 
gap in the minority spin channel. In this situation NQP 
states appear just above the Fermi level 

On the contrary, it was predicted that in half-metallic 
materials with a gap in the majority (say, "up") spin 



channel, NQP states should appear below the Fermi 
level 0, HI, E|- This asymmetry can be understood in 
terms of electron-magnon scattering processes, as pre- 
sented in the followings. 

A well studied model which takes into account the in- 
teraction of charge carriers with local moments is the s-d 
exchange model. The interacting part of the Hamilto- 
nian is given by —IJ2 ^i a a^i a c i^i where I is the s-d 
exchange parameter, Si represents the localized spin op- 
erators, a a p are the Pauli matrices, and Ci a are operators 
for conduction electrons. The NQP picture turns out to 
be essentially different for the two possible signs of the 
s — d exchange parameter. 

The ground state of the system with I > (assum- 
ing that the Fermi energy is smaller than the spin split- 
ting 2IS) has maximum spin projection and, thus, the 
minority-electron band should be empty. For this case 
(I > 0), NQP states in the minority spin gap develop 
as a superposition of the majority-electron states plus 
magnon states, and of the minority-electron states, so 
that the total spin projection of the system is conserved. 
As a result of this spin-polaronic effect, the minority- 
electron density of states has a tail corresponding to 
the virtual conduction electron spin-flip processes with 
magnon emission. However, these virtual flips are im- 
possible below the Fermi energy Ep due to the Pauli 
principle (all majority-electron states are already occu- 
pied and thus unavailable). Therefore, for the positive 
s-d exchange interaction the NQP states form above Ep. 

Contrary, for negative / the minority-spin band lies be- 
low the majority-spin one @, 0|. Occupied minority-spin 
states can be superposed with majority-electron states 
plus magnons, with conserved total spin projection, so 
the NQP states occur below Ep. At the same time, for 
I < the ferromagnetic ground state is non-saturated 
and thus zero-point magnon fluctuations are allowed. It 
is the fluctuations which are responsible for formation of 
occupied majority-electron NQP states there. 

Formally, the difference between I < and the previ- 
ous / > cases, can be explained in terms of a particle- 
hole transformation c\ a — > dig , and Cj CT — > d is . This mod- 
ifies the s-d exchange Hamiltonian into I ^2 SiO- a pd\ a dip. 



2 



In other words, the Hamiltonian with / > for electrons 
is equivalent to that with I < for holes. 

The above argument based on the s-d exchange model 
can be generalized for arbitrary multi-band half-metallic 
electronic structures 0, Q. The conclusion remain un- 
changed: for the case of minority- electron gap, NQP 
states are situated above the Fermi energy, while for 
the cases when the gap is present for majority-electrons, 
NQP states are formed below the Fermi energy. 

Most HMF materials have a gap in the minority spin 
channel so that NQP states arise above the Fermi level. 
As a consequence, these states cannot be studied by 
the very well-developed and accurate technique of spin- 
polarized photoemission [l9| , which can only probe occu- 
pied states. The spin-polarized Bremsstrahlung Isochro- 
mat Spectroscopy (BIS) probing unoccupied states [2(| 
has a much lower resolution. For this reason, HMF with 
a gap in the majority-spin channel, and, consequently, 
NQP states in the occupied region of the spectrum, allow 
for a detailed experimental analysis of these correlation- 
induced states and are, therefore, potentially of great in- 
terest. It is the purpose of the present work to perform an 
electronic structure calculation based on a combination 
of the generalized-gradient approximation (GGA) and of 
DMFT for the half-metallic fcrrimagnetic full-Heusler al- 
loy M112VAI, which has a gap in the majority spin chan- 
nel. By appropriately taking into account effects due to 
electronic correlations, we demonstrate explicitly the ex- 
istence of majority spin NQP states arising just below 
the Fermi level, and study the temperature dependence 
of their spectral weight. 

In full Hcusler compounds with the formula X2YZ, Mn 
atoms usually occupy the Y-position, while compounds in 
which Mn assumes the X-position Mn2YZ, are very rare. 
The prototype from the latter category is M112VAI, for 
which a large number of theoretical and experimental in- 
vestigations have been made. Neutron diffraction experi- 
ments [1H demonstrated the existence of a ferrimagnetic 
state in which Mn has a magnetic moment of 1.5±0.3/ib 
and V moment is — 0.9ug. X-ray diffraction and magneti- 
zation measurements j22[ found a total magnetic moment 
of 1.94/is at 5K, close to the half-metallic value of 2/i£. 
The Curie temperature of the sample was found to be 
about 760K and the loss of half-metallic character was 
attributed to the small amount of disorder. Electronic- 
structure calculations performed by Ishida [23| within the 
local-density approximation (LDA), predict the ground 
state of Mn2VAl to be close to half-metallicity. Weht 
and Pickett [HI used the GGA for the exchange cor- 
relation potential and showed that M112VAI is a half- 
metallic ferrimagnet with atomic moments in very good 
agreement with the experiment. Recent calculations of 



The intermixing between V and Al atoms in the M112VAI 
alloy showed that a small degree of disorder decreases the 
spin polarization at the Fermi level from its ideal 100% 
value, but the resulting alloy Mn^Vi^, Ali^ still show 
an almost half- metallic behavior [26|, |27| • 

According to the ideal full Heusler (L2i) structure, the 
V atom occupy the (0, 0, 0) position, the Mn atoms are 
situated at (1/4, 1/4, l/4)a and (3/4, 3/4, 3/4)a, and the 
Al at (1/2,1/2,1/2)0, where a = 5.875A is the lattice 
constant of the Mn2VAl compound. In our work, corre- 
lation effects in the valence V and Mn d orbitals are in- 
cluded via an on-site electron-electron interaction in the 
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the exchange parameters for M^VAl [25[ show a strong 
M n — V exchange interaction that influence the ordering 
in the Mn sublattice. The estimated Curie temperatures 
are in good agreement with the experimental values [25| . 



form I 

interaction is treated in the framework of dynamical 
mean field theory (DMFT) [l3j] , with a spin-polarized T- 
matrix Fluctuation Exchange (SPTF) type of impurity 
solver (28j. Here, Ci ma /c\ ma destroys/creates an elec- 
tron with spin a on orbital m on lattice site i. The 
Coulomb matrix elements U mm >m"m"> are expressed in 
the usual way [29] in terms of three Kanamori parame- 
ters U, U' = U — 2 J and J. Typical values for Coulomb 
(U = 2eV) and Stoner (J = 0.93eV) parameters were 
used for Mn and V atoms. The above value of U is con- 
siderably smaller than the bandwidth of Mn2 VA1 (7-8eV) 
therefore the use of a perturbative SPTF-solver is justi- 
fied. In addition, the same solver was used to investi- 
gate spectroscopic properties of transition metals with 
remarkable results 

Since the static contribution from correlations is al- 
ready included in the local spin-density approximation 
(LSDA/GGA), so-called "double counted" terms must 
be subtracted. To achieve this, we replace T, a (E) with 
Y, a (E) - £ CT (0) [HI in all equations of the DMFT pro- 
cedure [lj|. Physically, this is related to the fact 
that DMFT only adds dynamical correlations to the 
LSDA/GGA result. For this reason, it is believed that 
this kind of double-counting subtraction "£(0)" is more 
appropriate for a DMFT treatment of metals than the 
alternative static "Hartree-Fock" (HF) subtraction [36l |. 

In Fig. Q] we present the total density of states com- 
puted in GGA and GGA+DMFT, for T=200K. The 
GGA density of states displays a gap of about 0.4eV in 
the majority spin channel in agreement with previous 
calculations [24j. As expected from the s-d model calcu- 
lation, majority spin NQP states are visible just below 
the Fermi level, with a peak around — 0.25ey. In order 
to evaluate the spectral weight of these NQP states, we 
fit the low-energy density of states below the Fermi level 
in the majority channel with a Gaussian centered around 
the peak position. The NQP spectral weight is then de- 
fined as the area below the Gaussian curve. The inset 
shows the NQP spectral weight for several temperatures 
up to 300K. It is interesting to note that within the com- 
puted temperature range 50 < T < 300 these values are 
almost constant and considerably larger in comparison 
with similar values for (NiFe)MnSb [16[. The data pre- 
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FIG. 1: (color online) Total density of states, computed 
within GGA (dashed/blue), and GGA+DMFT (full/red). 
The Gaussian fit to the density of NQP states is shown as a 
dotted-dashed (black) line just below the Fermi level for the 
majority spin channel. The temperature dependent spectral 
weight of NQP states is displayed in the inset. 
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FIG. 2: (color online) Atom resolved density of states, com- 
puted within the GGA (dashed/blue) and GGA+DMFT 
(full/red) approach, at T=200K. The majority spin NQP 
states are visible in the Mn-3d T DOS just below the Fermi 
level. 



sented in the inset can be extrapolated down to T = OK, 
and a spectral weight of w 0.544±0.018(siaies/Afn— d) is 
obtained. This demonstrates that NQP states are present 
also at T=0K, and obviously they are not captured by the 
mean-field, GGA result. As it will be discussed below, 
NQP states predominantly consist of Mn-d electrons, so 
the value for the integrated spectral weight (inset of Fig. 
[I} only comes from Mn-d orbtials. 

The atom resolved DOS is presented in Fig. fJl In 
GGA, the net magnetic moment per unit cell is 2/i# with 
parallel Mn moments having values close to 1.6/ib and 
oppositely oriented V moments close to — 0.8^s. Be- 
low the gap, the majority spin total DOS is mainly of 
Mn character. The Mn and V moments have a strong 
t2 g character, and a small Al contribution to the mag- 
netic moment is present. Most of the V majority spin 
states lie above the gap, along with the Mn e g states. 
Minority-spin states below 0.5eV have roughly an equal 
amounts of V(t 2g ) and Mn character. States around the 
Fermi energy have a predominant Mnit^g) character, in 
agreement with [24]]. In contrast to the GGA results, 
the many-body DMFT calculation (see Fig. [2]) yields a 
significant DOS for the majority spin states just below 
the Fermi level. These are the majority spin NQP states 
discussed above [l2, 14, 15, 1(| 13, E3- As can be seen 
in Fig. [5] majority spin NQP states are predominantly 
of Mn — id) character. Their spectral weight is quite 



significant (see inset of Fig. [T]) so that accurate spin- 
polarized photoemission experiments should be able to 
identify the existence of such states. In contrast, ma- 
jority spin V(t2 S ) states below the Fermi energy are not 
significantly changed. Above the Fermi level, the Mn(e 9 ) 
and V(e 9 ) states are pushed to higher energy, such that a 
gap is formed just above Ep . In the minority spin chan- 
nel below Ep, both V and Mn(t2 S ) states are slightly 
modified, while above Ep, V(e s ) states are shifted to 
higher energies by O.beV. Around the Fermi level, the 
dominant Mn — 3d^ DOS is not significantly changed 
with respect to the GGA values. 

The applicability of the local DMFT approach to the 
problem of the existence of NQP states has been dis- 
cussed in ref. and 0. It is essential to stress that 
the accurate description of the magnon spectrum is not 
important for the existence of nonquasiparticle states and 
for the proper estimation of their spectral weight, but can 
be important to describe an explicit shape of the density 
of states "tail" in a very close vicinity of the Fermi energy. 

The imaginary part of the atom and orbital resolved 
self-energies, for T=200K are presented in Fig. [3] For the 
minority Mn-e g , V-e g and V-f2g-orbitals we observe that 
the imaginary part of the self-energy has a rather sym- 
metric energy dependence around the Fermi level, with 



a normal Fermi-liquid-type behavior -JmE 
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FIG. 3: (color online) Imaginary part of the self-energies 
lraYf Mn i V for t2 9 -orbitals on Mn (left upper panel) and V 
(right upper panel). The solid (red) line shows results for the 
majority(l) spins, while the dashed (blue) line for minority 
spins. The lower pannels show the corresponding 7m£^ /n / v 
for the e 3 orbitals. 



(E — Ep) 2 . The majority spin — ImY? M > v (E) shows 
a significant increase right below the Fermi level which 
is more pronounced for the £2g-orbitals. In addition, a 
slight kink is evidenced for an energy around -0.25eV. 
The majority-state nonquasiparticles are visible in the 
majority spin channel Fig. [5] at about the same energy. 
These results shown in Fig. [3] suggests that many-body 
effects are stronger on Mn than on V sites. Therefore, 
NQP states are mainly determined by the Mn-d atoms. 

The behavior of the imaginary part of the self-energy 
(Fig. [3]) and the Green function (Fig. [2]) can be corre- 
lated with the analysis of the spin-resolved optical con- 
ductivity. We have estimated the latter for different 
temperatures whithin the GGA and GGA+DMFT ap- 
proaches using an approximation of constant matrix ele- 
ments. Already at 50K the majority-spin optical spectra 
shows the appearance of a Drude peak signaling the clo- 
sure of the majority spin gap. With increasing temper- 
atures, spectral weight is transfered towards the Drude 
peak contributing to the depolarization discussed below 
in Fig. |U 

As it was demonstrated previously @, 0, HI, @ , the non- 
quasiparticle spectral weight in the density of states (Fig. 
[2]) is proportional to the imaginary part of the self-energy 
(Fig. [3|), therefore it is determined by the quasiparticle 
decay, which is the reason for the name of these states. 

Fig. 2] displays the temperature dependence of the 



FIG. 4: (color online) Temperature dependence of the spin 
polarization of conduction electrons (full/red) at the Fermi 
level P(Ef,T), and normalized magnetization M(T)/M(0) 
(dashed/black). 



magnetization obtained directly from the GGA+DMFT 
calculations and the spin polarization at the Fermi level, 
obtained using the relation P(E F ,T) = (N^(E F ,T) - 
N l {E F ,T))/(N ] (E Fl T)+N l {E F ,T)), N a being the den- 
sity of states. These results reflect a general trend 
valid for half-metals in the presence of NQP states 
0, ED, El El El EE El, namely that magnetization 
and polarization behave differently as function of temper- 
ature. However, as can be seen in Fig. 2J this difference 
is not as sharp as in other HMF materials [HE!- 

The effect of disorder on half-metallicity was recently 
discussed in Mn 2 Vi_a;Ali +a ; alloys, for —0.2 < x < 0.2 
H S E3- The excess of both Al and V atoms (x 
0. 1 - 0.1 or x = 0.2/ - 0.2) has the effect of shrinking 
the gap to zero, but with the Fermi level situated within 
the gap. In addition, the Mn moment is not affected by 
disorder and remains constant, in contrast to the V mo- 
ment. Spin polarization is decreased by about 10%, with 
electrons around the Fermi level having a dominant mi- 
nority spin character [26| . In contrast, many-body corre- 
lations have a much more dramatic effect. For non-zero 
temperatures, all atomic magnetizations are decreased. 
For instance, near room temperature (T = 300-K") the 
strongest decrease occurs in V, for which the moment 
drops almost by 47%, the Mn moment is reduced by 
32%, while the Al experiences just a small reduction by 
4%. As one can see from Fig. SI already at 5QK polar- 
ization drops to 75%, and is further decreased upon in- 
creasing the temperature. As we discussed previously for 
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the case of FeMnSb [16| , many-body induced depolariza- 
tion is significantly stronger than the effect of disorder or 
of other spin-mixing mechanisms such as spin-orbit cou- 
pling. This observation seems to hold also for the case of 
Mn2VAl, although we can not exclude the fact that for 
a larger degree of disorder, the material could possibly 
depart from its almost half-metallic situation displayed 
for small degree of substitution (—0.2 < x < 0.2). 

In conclusion, in this paper we have shown for a specific 
material that NQP states are also present in half-metals 
with a gap in the majority spin channel, and appear 
just below the Fermi level, as predicted in model cal- 
culations In the case of M112VAI, these states mainly 
consist of Mn-SS electrons and have a considerable spec- 
tral weight. Although this material was reported to be a 
half- metal from electronic structure calculations 24J , the 
experimental evidence is not clear. Several reasons are 
invoked such as existence of defects or the reduced sym- 
metry at surface and interfaces. From a theoretical point 
of view, we show that correlation- induced NQP states sig- 
nificantly change the majority spin electronic states, thus 
reducing the spin polarization at Ep. The appearance of 
NQP states and its connection with tunneling magne- 
toresistance was recentlystudied in Co2MnSi-based tun- 
nel magnetic junction [l2|. A great challenge would be 
to produce TMR junctions based on the ferri-magnetic 
Mn2VAl. This would allow for a direct experimental in- 
vestigation of the existence of majority spin NQP states. 
Promising candidate HMF materials with a majority spin 
gap of similar magnitude as M^VAl are the double per- 
ovskites Sr 2 FeM0 6 (M=Mo,Re) or Sr 2 CrRe0 6 often as- 
sociated with collosal magnetoresistance behavior. In 
particular the electronic structure of Sr 2 CrRe06 shows 
a closure of its maj ority spin gap in the presence of spin- 
orbit coupling [37], with states having a small spectral 
weight symmetrically distributed around the Fermi en- 
ergy. Discrepancies between the experiment and theoret- 
ical computations were explained based on possible anti- 
site disorder [37J. We suggest that a significant density 



of NQP states could be present in the above perovskites 
as well. Work on these lines is in progress. 

L.C. and E.A. acknowledge financial support by the 
Austrian science fund under project nr. FWF P18505- 
N16. L.C. also acknowledge the financial support offered 
by Romanian Grant CNCSIS/ID672/2009. M.I.K. ac- 
knowledges financial support from FOM (The Nether- 
lands). A.I.L. acknowledge financial support from the 
DFG (Grants No. SFB 668- A3). 



[1] 
[2] 



R. A. de Groot, F. M. Mueller, P. G. van Engen, and 
K. H. J. Buschow, Phys. Rev. Lett. 50, 2024 (1983). 
M. I. Katsnelson, V. Y. Irkhin, L. Chioncel, A. I. Licht- 
enstein, and R. A. de Groot, Reviews of Modern Physics 
80, 315 (2008). 



[3 
[4 
[5 
[6 
[7; 
[8 
[9 
[10; 
[11 
[12 

[13 

[14 

[is; 

[16 

[17; 

[is; 

[19 
[20 
[21 

[22 

[23 

[24 

[25 

[26 

[27 

[28 

[29 

[30 

[31 



I. Zutic, J. Fabian, and S. D. Sarma, Rev. Mod. Phys. 
76, 323 (2004). 

D. M. Edwards and J. A. Hertz, Journal of Physics F- 
Metal Physics 3, 2191 (1973). 

V. Y. Irkhin and M. I. Katsnelson, J. Phys.: Condens. 
Matter 2, 7151 (1990). 

V. Y. Irkhin and M. I. Katsnelson, Phys. Usp. 37, 659 
(1994). 

V. Y. Irkhin and M. I. Katsnelson, Eur. Phys. J. B 30, 
481 (2002). 

E. McCann and V. I. Fal'ko, Phys. Rev. B 68, 172404 
(2003). 

G. Tkachov, E. McCann, and V. I. Fal'ko, Phys. Rev. B 
65, 024519 (2001). 

E. McCann and V. I. Fal'ko, Phys. Rev. B 66, 134424 
(2002). 

V. Y. Irkhin and M. I. Katsnelson, Phys. Rev. B 73, 
104429 (2006). 

L. Chioncel, Y. Sakuraba, E. Arrigoni, M. I. Katsnelson, 
M. Oogane, Y. Ando, T. Miyazaki, E. Burzo, and A. I. 
Lichtenstein, Phys. Rev. Lett. 100, 086402 (2008). 

G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, 
O. Parcollet, and C. A. Marianetti, Rev. Mod. Phys. 78, 
865 (2006). 

L. Chioncel, M. I. Katsnelson, R. A. de Groot, and A. I. 

Lichtenstein, Phys. Rev. B 68, 144425 (2003). 

L. Chioncel, M. I. Katsnelson, G. A. de Wijs, R. A. 

de Groot, and A. I. Lichtenstein, Phys. Rev. B 71, 085111 

(2005). 

L. Chioncel, E. Arrigoni, M. I. Katsnelson, and A. I. 
Lichtenstein, Phys. Rev. Lett. 96, 137203 (2006). 
L. Chioncel, P. Mavropoulos, M. Lezaic, S. Blugel, E. Ar- 
rigoni, M. I. Katsnelson, and A. I. Lichtenstein, Phys. 
Rev. Lett. 96, 197203 (2006). 

L. Chioncel, H. Allmaier, A. Yamasaki, M. Daghofer, 
E. Arrigoni, M. Katsnelson, and A. Lichtenstein, Phys. 
Rev. B 75, 140406 (2007). 

P. D. Johnson, Rep. Prog. Phys. 60, 1217 (1997). 
M. Donath, Surf. Sci. Rep. 20, 251 (1994). 

H. Itoh, T. Nakamichi, Y. Yamagichi, and N. Kazama, 
Trans. Japan. Inst. Met. 24, 256 (1983). 

C. Jiang, M. Venkatesan, and J. M. D. Coey, Solid State 
Commun. 118, 513 (2001). 

S. Ishida, S. Asano, and J. Ishida, J. Phys. Soc. Jpn. 53, 
2718 (1984). 

R. Weht and W. E. Pickett, Phys. Rev. B 60, 13006 
(1999). 

E. Sasioglu, L. M. Sandratskii, and P. Bruno, J. Phys.: 
Condens. Matter 17, 995 (2005). 

K. Ozdogan, I. Galanakis, E. Sasioglu, and B. Aktas, J. 
Phys.: Condens. Matter 18, 2905 (2006). 

I. Galanakis, K. Ozdogan, E. §a§ioglu, and B. Aktas, 
Phys. Rev. B 75, 092407 (2007). 

M. I. Katsnelson and A. I. Lichtenstein, Eur. Phys. J. B 
30, 9 (2002). 

M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

J. Minar, L. Chioncel, A. Perlov, H. Ebert, M. I. Kat- 
snelson, and A. I. Lichtenstein, Phys. Rev. B 72, 045125 
(2005). 

J. Minar, H. Ebert, C. D. Nadai, N. B. Brookes, F. Ven- 
turing G. Ghiringhelli, L. Chioncel, M. I. Katsnelson, and 
A. I. Lichtenstein, Phys. Rev. Lett. 95, 166401 (2005). 



6 



[32] J. Braun, J. Minar, H. Ebert, M. I. Katsnelson, and A. I. 
Lichtenstein, Phys. Rev. Lett. 97, 227601 (2006). 

[33] S. Chadov, J. Minar, H. Ebert, A. Perlov, L. Chioncel, 
M. I. Katsnelson, and A. I. Lichtenstein, Phys. Rev. B 
74, 140411 (2006). 

[34] A. Grechnev, I. Di Marco, M. I. Katsnelson, A. I. Lichten- 
stein, J. Wills, and O. Eriksson, Phys. Rev. B 76, 035107 
(2007). 



[35] A. I. Lichtenstein, M. I. Katsnelson, and G. Kotliar, 
Phys. Rev. Lett. 87, 067205 (2001). 

[36] A. G. Petukhov, I. I. Mazin, L. Chioncel, and A. I. Licht- 
enstein, Phys. Rev. B 67, 153106 (2003). 

[37] G. Vaitheeswaran, V. Kanchana, and A. Delin, Appl. 
Phys. Lett. 86, 032513 (2005). 



